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Scattering rate calculations in two-dimensional Si/Sii-xGca; systems have typically been restricted 
to rectangular Ge profiles at interfaces between layers. Real interfaces however, may exhibit diffuse 
Ge profiles either by design or as a limitation of the growth process. It is shown here that alloy 
disorder scattering dramatically increases with Ge interdiffusion in (100) and (111) n-type quantum 
wells, but remains almost constant in (100) p-type heterostructures. It is also shown that smoothing 
of the confining potential leads to large changes in subband energies and scattering rates and a 
method is presented for calculating growth process tolerances. 

PACS numbers: 73.50.Bk, 73.21.Fg, 73.63.Hs, 73.61. Cw 

Keywords: Silicon; germanium; SiGe; carrier dynamics; intersubband transitions; interface roughness; alloy 
disorder; Coulombic interactions; phonon interactions 



I 

O 

o 



> 
p 

00 

o 

o 

> 

• I— I 

X 
S3 



I. INTRODUCTION 

Two dimensional intersubband devices in the 
Si/Sii_a;Gea: materials system offer a possible reduction 
in fabrication costs compared with more conventional 
III-V systems.^ Successful operation of resonant tun- 
neling diodes (RTDs) has been achievable for several 
years and electroluminescence from more complex 
quantum cascade structures has been observed^''-- A 
quantum cascade laser (QCL) has not yet been devel- 
oped in Si/Sii_2;Ge2; although several designs have been 
proposed recentlyi^i^ 

In order to accurately design and simulate such struc- 
tures, a good understanding of intersubband carrier dy- 
namics is required. Previous models have assumed that 
interfaces between layers are perfectly abrupt?^ while in 
reality, diffuse Ge profiles may result either by design, 
by interdiffusion during the growth process or by surface 
segregation of Ge atoms4 

The Ge interdiffusion dramatically changes the sub- 
band spacing^ and the overlap between wave functions. 
A more accurate model of intersubband scattering rates 
must therefore account for these effects. In this paper, 
we review the models of the principal intersubband scat- 
tering rates and extend the conventional interface rough- 
ness scattering model to an arbitrary interface geometry. 
To determine the effect on simple intersubband systems, 
scattering rates were calculated as a function of subband 
spacing, electron temperature, and diffusion length in 
single quantum wells (QWs) in (100) p-type and (100) 
and (111) n-type systems. 

Knowing the effect of interdiffusion on scattering rate, 
it is possible to estimate the robustness of device de- 
signs. By calculating intersubband scattering rates and 
subband spacings, design tolerance to interdiffusion may 
be estimated. The viability of a design may be assessed 
by comparing this with the capabilities of growth pro- 
cesses. We use the example of a coupled QW system to 
demonstrate the technique and predict the level of inter- 



diffusion required to cause device failure. 



II. SCATTERING MODELS 

Carrier scattering in a two dimensional Si/Sii_a:Gea; 
system may be described by a set of independent pro- 
cesses. The models for the Coulombic, electron-phonon, 
and alloy disorder interactions apply to arbitrary inter- 
face geometries, while some modification is required for 
interface roughness scattering. The mechanisms are sum- 
marized as follows. 



A. Interface roughness scattering 

In z-confined two-dimensional heterostructures, the 
confining potential varies with fluctuations in interface 
location over the [x, y) planej ^°i^^ The roughness is usu- 
ally assumed to have a Gaussian Fourier transform ^z{f) 
with height A and correlation length A, which is isotropic 
across the {x, y) plane jiSiii^ii* such that 



(A.(r)A,(r')) 



A2 



exp 



,/|2 



A2 



(1) 



The commonly used expression for the resulting scat- 
tering rate^^ assumes an abrupt interface geometry. This 
has been accurately fitted to experimental data for struc- 
tures with approximately abrupt interfacesj^^ but the ex- 
pression is incompatible with smooth envelope potentials. 
We determine the scattering rate for an arbitrary inter- 
face geometry and verify that it reduces to the specific 
case of an abrupt interface. 

The perturbation Av{R) due to a position shift A^(r) 
in an arbitrary confining potential V{z) is assumed to be 
correlated over the length of a single interface. At the 
point r = X + y assuming isotropy across the xy plane. 



Av{R) = V[z-AAr)]~V{z) 
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If the /-th interface in a multilayer structure is cen- 
tered about the plane z — zj and extends over the range 
(zlj, zuj), we define the scattering matrix element as 
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where |/) and \i) are the final and initial wave functions 
respectively and the rectangular function rect(z) is de- 
fined as 



-^(^) = { 0; \z\ ; oi • 
From Fermi's golden rule, the scattering rate Wfij{ki) 
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where kfi — kf — ki is the scattering vector, kf and ki 
are the final and initial wave vectors respectively, 9 is the 
scattering angle and md is the density-of-states effective 
mass. 

At this stage, the general result may be tested against 
the example of an abrupt change in envelope potential of 
magnitude Vb,/, where the perturbating potential is 



Av{R) = Vo(z/)rect 
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Under time-independent perturbation theory, the pertur- 
bation must be small, i.e. Az(r) and the perturbing 
potential simplifies to 

Av{R) = Vo{zi)A,ir)Siz-zi). (7) 

The scattering matrix element (Eqn. [3]) becomes 

Ff,j = {f\Vo,iS{z-zi)\i), (8) 

in agreement with the well known expressions^ 

Returning now to the general expression for scattering 
rate (Eqn.O, we determine the total scattering rate for a 
structure with N layers, numbered / € N"*",/ < N. The 
integral over 9 simplifies to a regular modified cylindrical 
Bessel function of zeroth order, Iqi^ By assuming rough- 
ness profiles are uncorrelated over separate interfaces, the 
total rate is found as a summation. 



e(a2). 
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where the Heaviside step function, O permits a nonzero 
rate only for real final wave vectors. The new matrix 
element is 



Bfi{ki 



where a 
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The new model was fitted to the experimental data de- 
scribed in detail in our previous paper »^ The parameters, 
{A=1.4A, A=50 A} accurately fit the measurements and 
are very similar to our previous theoretical values, and 
to other recent dataJ*^!^ The slight difference in fitting 
parameters arises from numerical approximations in the 
perturbing potentials: the perturbation for arbitrary in- 
terface geometries (Eqn. ^ uses a Taylor series expan- 
sion, whereas the solution for abrupt interfaces (Eqn. [T]) 
uses the Dirac S function limit of the narrow rectangular 
function. 



B. Alloy disorder scattering 



A standard scattering model for alloy disorder 



scattering 



19,20 



has been modified slightly to permit vari- 



able alloy composition, x{z). The resulting scattering 
potential becomes 



^}iz)x{z)[l-x{zM{z)dz, 

(11) 

where ao is the in-plane lattice constant and dV is com- 
monly approximated as the difference in conduction band 
potentials between Si and Ge. As a;(l — a;) = for a pure 
Si layer, the integral domain in abrupt interface systems 
is restricted to the barriers in n-type systems or the wells 
in p-type systems. In a diffuse system however, the Ge 
content is always nonzero and the integral domain ex- 
tends over the entire structure. 



C. Carrier— phonon scattering 

Wells in p-type heterostructures contain similar frac- 
tions of Si and Ge. Hole-phonon scattering via the de- 
formation potential interaction was therefore calculated 
for the Si-Si, Si-Ge, and Ge-Ge branches of the nonpo- 
lar optical modei^i^i In n-type heterostructures however, 
the Ge fraction in the wells is small and only the Si-Si 
branch was considered. Intravalley deformation potential 
scattering for the acoustic mode was included for both n- 
and p-type heterostructures. 

In n-type, Si-rich systems, scattering between conduc- 
tion band A valleys is described by either g processes, 
which transfer electrons to the opposite valley in recip- 
rocal space or / processes which transfer electrons to the 
four perpendicular valleysi^ ff-LO, /-LA, and /-TO in- 
teractions are permitted in a zero-order model, whereas 
f/-LA, f;-TA, and /-TA interactions are permitted only 
as first-order processes^ and are somewhat slower. The 
interactions are characterized by a deformation potential 
and a frequency, loq, which for intervalley interactions is 
nonzero. The scattering rates increase rapidly with sub- 
band separation, Efi until they saturate at > Hloq. 
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D. Coulombic interactions 

Ionized impurity scattering and carrier-carrier scat- 
tering rates were calculated as Coulombic interactions 
between either a carrier and a dopant ion or a pair 
of carriers. The expression for ionized impurity scat- 
tering given by Unumafi^ was modified to incorporate 
static screening in the Thomas-Fermi approximation^^'* 
Carrier-carrier scattering rates were calculated using the 
screened Coulombic interaction model described by Smet 
et. al.^^ Both rates are fastest for transitions between en- 
ergetically similar states. Doping was set as 1 x 10*® cm^"^ 
throughout each structure we considered, as modulation 
doping is difficult to achieve in SiGe heterostructures4 
We also assume that all dopants are ionized at low tem- 
peratures. The sheet doping density was calculated as 
the product of the total volume doping and the length of 
the structure. 



E. Average rates 

As justified previouslyj-- all subband electron temper- 
atures were set to a single average value, Tg, assumed 
to be different from the lattice temperature, which was 
taken as T—AK in our calculations. The average scat- 
tering rate from the second to first subband, W12 was 
calculated as 



vyi2 = , 



(12) 



where intrasubband scattering was assumed to be much 
faster than intersubband. The initial distribution of elec- 
trons is therefore an equilibrium Fermi-Dirac function, 
fl^i^^) using the quasi- Fermi level for the subband, 
where fc2 is the initial wave vector. The assumption has 
also been made that the destination states are always 
unoccupied, which is reasonable at low doping levels. 



III. DIFFUSE QUANTUM WELLS 

Annealing of an abrupt structure, with the Ge frac- 
tion xi in layer / provides a simple model of a diffuse 
system. The abrupt-interface structure is embedded be- 
tween infinitely thick barriers with composition a;o- The 
composition profile after annealing is^ 



1 ^ 



XI 



7=1 



erf 



Zl-l 



Xo^ 

2 



erf 



erf 



erf 



z- zo 



z — Zj 

L 



(13) 



where the /-th layer boundaries are zj-i and z/, and L 
is the diffusion length. 

In the calculations presented below, we have assumed 
that the composition profiles are symmetrical for the left 
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FIG. 1; A2 conduction band edge for a 10 nm Si QW be- 
tween two 5nm Sio.sGeo.s barriers, with varying Ge diffusion 
lengths. The A4 conduction band edge for an abrupt interface 
is shown to be 120 meV higher in energy. 



and right interfaces of a QW. This corresponds to the 
case where interdiffusion dominates over surface segre- 
gation. If Ge segregation effects are strong, then the 
interface profiles will be asymmetric, although the effect 
on scattering is expected to be similar. 

A 10 nm (100) n-type QW between two 5 nm Sio.sGeo.s 
barriers is shown in fig. [H with the in-plane lattice con- 
stant set to achieve strain balance. As L increases, the 
bottom of the well narrows and the top widens. The ef- 
fect on scattering in (100) n- and p-type and in (111) 
n-type systems is discussed in depth in the following sec- 
tions. 

Interdiffusion increases the separation of low energy 
subbands as shown in Fig. [51 As scattering rates depend 
on subband separation, interdiffusion affects scattering 
in two possible ways: 

1. A direct effect due to the change in interface geom- 
etry. 

2. An indirect effect due to the change in subband 
spacing. 

Throughout this section, we adjust the width of QWs to 
correct the interdiffusion effect on subband separation. 
The calculated change in scattering rates is therefore due 
only to the change in interface geometry. In section ITVl 
the total effect is determined by varying interdiffusion 
without correcting the subband separations. 



A. (100) n-type single QW 

In Sii_a;Gea; alloys with x < 85%, the conduction band 
has six minima near the Brillouin zone edge in the A 
directions»2I The valleys are almost parabolic, and a sin- 
gle band effective mass approximation (EMA) accurately 
models electron confinement 



4 




Well width (nm) Electron temperature (K) 



FIG. 2: Energy separation between the lowest pair of sub- 
bands in a Si/Sio.sGeo.s QW as a function of well width. Re- 
sults are shown for A valleys in n-type structures (Si wells) in 
the (100) and (111) orientations and for p-type structures (Si 
barriers) in the (100) orientation. In all cases, the results are 
shown for interface diffusion lengths of Onm, 1 nm and 2nm 
with arrows denoting the increasing diffusion length. 



The A valleys have ellipsoidal equipotential surfaces 
with their major axes along the cubic unit cell edge di- 
rections. Two separate effective masses are defined for 
electrons with wave vectors near the A minima: the lon- 
gitudinal mass, mi = 0.916 and the transverse mass, 
rrit = 0.19.— Quantum confinement in the z direction 
in strained (100) systems yields fom degenerate A4 sub- 
bands with quantization effective mass, rrig = rrif and two 
A2 subbands with rriq — mi. As shown in fig. [l] the A4 
conduction band edge is at a relatively high energy and it 
can be assumed to have a negligible electron population. 
The following discussion therefore considers only scatter- 
ing between A2 subbands and omits /-phonon emission 
processes. 

A more precise double valley EMA models states with 
mq = mi as a combination of basis states from the two 
A2 valleys. In z confined systems, the phase difference 
between reflected basis components splits the degeneracy 
of the subbands. We have previously shown however, 
that the splitting is quite small for QWs wider than 2- 
3nm and it is therefore omitted in this work.^'' 

The width of a Si QW between a pair of Sio. 5600.5 
barriers (as in fig. [T]) was adjusted for £'2i = 10meV for a 
given diffusion length. Scattering rates were then calcu- 
lated as a function of and are shown in fig. [3l 

All mechanisms except alloy disorder scattering were 
found to be almost independent of diffusion length and 
their rates are only plotted for the abrupt interface sys- 
tem, for simplicity. For all diffusion lengths, Coulom- 
bic interactions and interface roughness rates were rela- 
tively large. As intcrdiffusion increases however, the al- 
loy disorder rate increases very rapidly, from 5xl0^s~^ 
for L=Onm to LlxlG^^s"^ for L=2nm. As described 



FIG. 3: Average scattering rates from the second to first 
subband in a n-type (100) 10 nm Si QW between two 5nm 
Sio.sGeo.s barriers with r=4K as a function of electron tem- 
perature. Rates are shown for electron-electron (E_E), in- 
travalley acoustic phonon (AC), interface roughness (IFR), 
ionized impurity (ION), g-type phonon emission, and alloy 
disorder (ADO) processes. ADO is strongly dependent on 
diffusion length and is shown at L^Onm (no symbols), Inm 
(circles), and 2nm (stars). All other rates are almost inde- 
pendent of interdiffusion and are shown only for L=0 nm 



previously, this is due to alloy disorder scattering being 
permitted in the well region of diffuse structures, where 
the electron probability is large. 

The electron temperature dependence of the average 
alloy disorder and intravalley acoustic phonon scatter- 
ing rate is shown to be weak in fig. [3l as the scatter- 
ing potentials depend neither on the initial wave vector 
nor explicitly on the electron temperature. The inter- 
face roughness scattering however depends on the initial 
wave vector and hence the average rate is affected by the 
temperature dependent distribution of electrons in the 
subband. The screening of the Coulombic interactions 
is also affected by electron temperature, explaining the 
gradual decrease in average rate. Finally, the intervalley 
phonon emissions are affected very strongly by electron 
temperature because the subband spacing is lower than 
huio for each of the permitted phononsj^ The average 
scattering rate therefore depends on the number of elec- 
trons with sufficient initial kinetic energy to scatter into 
a state within the lower subband. 

In pump-probe experiments, which are often used to 
determine scattering lifetimes, the carrier temperature, 
Te is elevated above the lattice temperature, T and de- 
cays over time. We have previously shown however, that 
good agreement with experimental data is achievable by 
assuming = T+20K.^- In fig. HI we show the average 
intersubband scattering rate as a function of subband 
separation at T =4K, re=24K. 

For L=Onm and subband spacing closer than lOmeV, 
ionized impurity and electron-electron scattering domi- 
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FIG. 4: Average scattering rates from second to first subband 
in a (100) n-type Si QW as a function of subband separation 
at electron temperature of 213=24 K. The ADO rate is shown 
at diffusion lengths, L=Onm (no symbols), Inm (circles), and 
2 nm (stars). All other rates are almost independent of L and 
are shown only for L=2 nm for simplicity. 



nate, while at spacings between 10 meV and 55 meV, in- 
terface roughness and intravalley acoustic phonon scat- 
tering are fastest. As the subband spacing becomes com- 
parable to the energy of the g-hO phonons, the emission 
rate exceeds the intravalley acoustic phonon scattering 
rate. Alloy disorder scattering is again shown to increase 
significantly with interdiffusion and becomes dominant 
for subband spacing above lOmeV and L=2nm. 



B. (Ill) n-type single QW 

Several important changes are introduced by moving 
to the (111) orientation. Electron confinement is now de- 
termined using an oblique cross-section of the ellipsoidal 
equipotential surface in fc-space, which is identical in all 
six A valleys. The quantization effective mass is^ 



Smirrit 
2m; + nit 



0.26, 



(14) 



compared with ~ m;=0.916 for the A2 valleys in 
(100) systems. The density-of-states effective mass iaS^ 



mt{2mi + mt) 



= 0.36, 



(15) 



compared with md = mt=G.19 for the A2 valleys in (100) 
systems. 

The uniaxial strain splitting exhibited in (100) is ab- 
sent in (111) systems, while the hydrostatic strain in- 
duced shift in conduction band potential is given by 
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FIG. 5: Average scattering rates from second to first subband 
in a (111) n-type Si QW as a function of subband separation 
with Te=24K. All rates from fig.|4]are shown (using the same 
legend) as well as the /-phonon rates. 



where cu are elastic constants, :iu,d are deformation po- 
tentials, and £11 is the in-plane strain as opposed to 



2i?il ( + gS„ 



cii - C12 
Cu 
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for (100) systems.— For a strain symmetrized 10 nm Si 
QW between two 5 nm Sio.sGeo.s barriers, the conduction 
band offset is around 150 meV, compared with around 
380 meV for the A2 offset in (100) as shown in fig. [T] As 
the valley degeneracy is not split, the population of sub- 
bands in each valley is equal and /-phonon processes are 
no longer negligible. Fig. [5] shows the average intersub- 
band scattering rates in a (111) n-type Si QW between 
two 5nm Sio.sGco.s barriers with Te=24K as a function 
of subband separation. 

The Coulombic interactions and the intravalley- and 
f/-phonon interactions are almost unchanged compared 
with the (100) case. From eqn. [51 it can be deduced that 
the reduced conduction band offset reduces the interface 
roughness scattering rate in (111) systems slightly. More 
significantly, a larger well width is required to achieve an 
equivalent subband splitting and hence a smaller propor- 
tion of the wave functions extend over the interface re- 
gion. The scattering matrix element is therefore smaller 
than in the (100) case. The alloy disorder scattering is 
also slightly increased by the change in dcnsity-of-states 
effective mass. 

The zero-order /-phonon emission processes become 
large at E12 > hwo and in structures with abrupt in- 
terfaces represent the dominant mechanisms, /-phonon 
interactions have four destination valleys and therefore 
dominate over (/-phonons, which have only a single des- 
tination. 
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FIG. 6: Average scattering rate from second to first HH sub- 
band in a p-type QW as a function of liole temperature, witii 
subband separation fixed at lOmeV. The legend is the same 
as previous figures, with the optical phonon branches as la- 
beled. The alloy disorder scattering rate is shown at diffusion 
lengths, L—Q nm (no symbols), L—1 nm (circles), and L=2 nm 
(stars). All other rates are shown at L=Onm only. 



FIG. 7: Average scattering rates from second to first subband 
in a (100) n-type double QW with Si wells of width 5 nm and 
3nm separated by a Inm Sio.sGeo.s barrier as a function of 
Ge interdiffusion. The legend is the same as that in fig.|4] The 
subband separation as a function of interdiffusion is overlaid 
(crosses) . 



C. (100) p-type single QW 

The band structure for hole transitions is highly non- 
parabolic and state contributions from the light hole 
(LH), heavy hole (HH), and spin split-ofF (SO) bands 
are all significant at nonzero in-plane wave vectors. The 
single band EMA is therefore inadequate and the 6x6 
k ■ p solution described previously)^ was used to account 
for the multiband effects. 

As Si/SiGe forms type II hcterointerfaces, the well and 
barrier compositions of the structure used for n-type sys- 
tems were reversed. The resulting HH band offset was 
350 meV. The well width was adjusted to give a lOmeV 
separation between the two lowest HH states at fcy = 0. 
Fig. [6] shows the average intersubband scattering rates 
for the system as a function of hole temperature, T/j. 

The results differ considerably from those of the n-type 
systems considered previously. As Th increases, the hole 
distribution spreads over a larger range of in-plane wave 
vector and HH states acquire a larger LH contribution. 
This affects both the scattering matrix element and the 
effective density of states, resulting in an increase in all 
rates. 

The carrier-carrier scattering rate depends on the over- 
lap between initial and final states for a pair of carriers 
and is strongly dependent on the in-plane wave vector of 
the involved states and hence on hole distribution over 
fc||. The Th dependence is therefore extremely strong 
and carrier carrier scattering dominates above Th=65K 
in this system. Below this temperature, alloy disorder 
scattering is dominant. 

In contrast with n-type systems, alloy disorder scat- 
tering in p-type systems is very weakly dependent on 



interdiffusion. This is because the magnitudes of wave 
functions are largest in the center of the QW, where the 
Ge fraction is barely affected by the interdiffusion. 

IV. GROWTH PROCESS TOLERANCE 

We have so far considered the effect of interdiffusion 
on intersubband transitions of known energies. Interdif- 
fusion due to growth processes however changes the sub- 
band separation as well as the scattering rates. There are 
therefore important implications for fabrication of inter- 
subband devices such as RTDs and QCLs. The following 
method uses both these effects to estimate the tolerance 
of a device design to undesired interdiffusion. 

A (100) n-type double QW system with two coupled 
Si wells of width 5 nm and 3 nm, separated by a 1 nm 
Sio.5Gco.5 barrier provides a relatively simple example 
similar to that used in our previous pump-probe inves- 
tigation of intersubband transition lifetimes in p-type 
materials;^ The system was surrounded by a pair of 
5nm thick Sio.5Geo.5 barriers as before. For simplicity, 
the symmetric approximation for interdiffusion was pre- 
served and the temperatures were fixed at T—4K and 
Te=24K. The scattering rates and subband separation 
are plotted in fig. [T] 

In the nominal structure, the layers are defined pre- 
cisely, restricting alloy disorder scattering to the barrier 
regions. The separating barrier is thick enough for the 
coupling to be very weak between states. The subband 
separation of 17.5 meV is too low for phonon emission 
to be significant, although large enough (and the ma- 
trix element small enough) for Coulombic interactions to 
be negligible. The dominant rate is therefore interface 
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FIG. 8: A2 conduction band edge in a two-well heterostruc- 
ture, nominally comprising 50% Ge barriers and Si wells. The 
left and right well widths are 5 nm and 3 nm respectively, and 
the separating barrier width is 1 nm. The potential is shown 
for diffusion lengths of Onm, 0.7 nm and 1.5 nm, which co- 
incide with the uncoupled, weakly coupled and single well 
operating regimes respectively. 

roughness scattering, which is also slowed by the small 
overlap between states. As the interdiffusion increases, 
the device behavior varies and can be characterised by 
the following operating regimes, which are illustrated in 
fig.l 

1. Uncoupled wells (L < 0.7 nm): When interdif- 
fusion is small, extra Ge in the wells increases al- 
loy disorder scattering. The overlap of states is 
decreased as the bottom of the barrier becomes 
thicker. The other scattering rates therefore de- 
crease slightly. Conversely, the bottoms of the wells 
narrow and subband spacing is increased. 

2. Weakly coupled wells (0.7 ^ L < 1.5 nm): 

For moderate interdiffusion, the barrier potential 
is substantially reduced, while the subband spacing 
is increased further. As the second subband min- 
imum approaches the barrier potential, the over- 
lap between states increases rapidly. All scattering 
rates therefore begin to rise extremely rapidly. 

3. Single well (L > 1.5 nm): For large interdiffu- 
sion, the second subband energy exceeds the barrier 
potential and the character of the system changes 
from that of a double QW to a single QW with a 
small central perturbation. The perturbation di- 
minishes as interdiffusion increases and the sub- 
band separation decreases towards the value for a 
wide, single well. The effect on scattering is slightly 
more complex than the previous regimes. 

The first-order electron-phonon rates remain al- 
most constant as the effect of the increased over- 
lap is countered by the decrease in subband sep- 



aration. The 5-LO phonon emission however, de- 
creases sharply due to its high frequency and zero- 
order character. Coulombic interactions continue 
to increase as the overlap between states increases 
and the subband separation decreases. Interface 
roughness scattering slows as the barrier potential 
diminishes and alloy disorder scattering continues 
to increase as the Ge content in the single well rises. 

The total rate in fig. [7] remains almost constant at dif- 
fusion lengths below 1 nm, although the subband spac- 
ing varies considerably above 0.25 nm. A multiple QW 
device such as a QCL might therefore operate success- 
fully with 1 nm interdiffusion if changes in transition en- 
ergy are acceptable. Above this tolerance however, there 
are catastrophic changes in the device and both transi- 
tion energies and lifetimes will be severely affected. This 
method may readily be applied to a wide variety of in- 
tersubband device designs as an estimator of robustness. 

V. CONCLUSION 

We have simulated annealed interfaces to model real 
systems more accurately than the rectangular well ap- 
proximation. The study of carrier dynamics has thus 
been extended to systems with diffuse interfaces (either 
by design or via the growth process) . Intersubband scat- 
tering rates have been investigated as functions of sub- 
band separation and electron temperature in (100) p-type 
and (100) and (111) n-type diffuse QWs. 

We have shown that for any given subband separa- 
tion and electron temperature, the alloy disorder scat- 
tering increases rapidly with diffusion length in n-type 
systems as Ge content in the well region increases. The 
alloy disorder scattering is a relatively minor effect in sys- 
tems with abrupt interfaces, while in diffuse systems, it 
can become the dominant mechanism. In p-type systems 
however, the effect is negligible as the Ge content in the 
center of the well remains large. 

In the case of Ge interdiffusion introduced during 
growth, this implies that very simple p-type system de- 
signs may be adjusted to preserve subband separations 
and carrier dynamic behavior if the achievable diffusion 
length is known. In n-type systems however, the carrier 
dynamics may be strongly affected even when the design 
is adjusted to preserve subband separations. 

The robustness of an example double (100) n-type QW 
design was examined by considering the combined ef- 
fect of Ge interdiffusion on both the subband separation 
and scattering rates. It was shown that a tolerance of 
L <0.25nm preserved the subband separation, while a 
less restrictive tolerance of L < 1 nm preserved the scat- 
tering rates. 

As group IV heterostructure epitaxy is less developed 
than that of III-V systems, we propose that future de- 
signs of complex devices such as QCLs in the Si/SiGe 
materials system should be tested for their robustness 
using this method before attempting growth. 
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